EUROPHYSICS LETTERS 
Europhys. Lett., (), pp. () 



Self-assembly in mixtures of amphiphilic polymers and 
surfactants 

H. Diamant and D. Andelman 
School of Physics and Astronomy 

Raymond and Beverly Sackler Faculty of Exact Sciences 
Tel Aviv University, 69978 Tel Aviv, Israel 

(received ; accepted ) 

PACS. 61.25Hq - Macromolecular and polymer solutions, polymer melts, swelling. 
PACS. 61.41+e - Polymers, elastomers, and plastics. 

PACS. 87.15Da - Physical chemistry of solutions of biomolecules, condensed states. 



Abstract. - We present a model for the joint self-assembly of amphiphilic polymers and small 
amphiphilic molecules (surfactants) in a dilute aqueous solution. The polymer is assumed to 
consist of a hydrophilic backbone and a large number of hydrophobic side groups. Preference 
of the surfactant to bind to hydrophobic microdomains along the polymer induces an effective 
attraction between bound surfactants. This leads to two distinct binding regimes depending on 
a single physical parameter, e, which represents the ratio between surfactant-polymer affinity 
and polymer hydrophobicity. For small e the binding is non-cooperative, whereas for large e it 
becomes strongly cooperative at a well-defined critical aggregation concentration. Our findings 
are in accord with observations on diverse experimental systems. 



A considerable experimental and theoretical effort has been devoted in recent years to the 
study of mixtures of polymers and amphiphilic molecules (surfactants) in solution jjj. Of 
particular interest and importance are systems where the polymer itself has amphiphilic 
qualities as they exhibit rich self-assembly behaviour [0. These systems offer, on one hand, 
numerous industrial applications, since their properties (e.g., rheology) can be conveniently 
tuned and their solubility in water makes them "environmentally friendly" . On the other 
hand, their complex amphiphilic nature poses a challenge to current theoretical investigations 
far beyond more conventional polymers in solution. Besides industrial applications, biological 
macromolecules, such as proteins, RNA and single strands of DNA, may be viewed as special, 
more complex cases of amphiphilic polymers. 

In the current work we focus on the case where a large fraction of hydrophobic side groups 
are chemically attached to a hydrophilic polymer backbone (fig. la). The properties of such 
macromolecules and their interaction with free surfactants in solution have been investigated 
in a series of experiments during the past years, revealing rich self-assembly behaviour |3[ ||. 
One of the notable findings is that the cooperativity of the surfactant-polymer binding, which 
is related to interactions between bound surfactants, sensitively depends on the polymer prop- 
erties (e.g., its hydrophobic side chains and the degree of ionisation of its backbone). Previous 
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Fig. 1. - (a) Schematic view of a section of a side-chain amphiphilic polymer with bound 
surfactants, (b) Schematic view of the system as described by our model. 

theories treating polymer-surfactant mixtures |5|, [}| f?J did not address the diverse behaviour 
of those experimental systems. Our aim in this Letter is to construct a simple model which 
yields a qualitative understanding of the thermodynamics of such systems while accounting, 
in particular, for the polymer-induced interactions between bound surfactants. We neglect, 
therefore, structural details such as the size of the hydrophobic side chains, steric effects, chain 
rigidity etc., which are microscopically important but not expected to qualitatively affect the 
thermodynamics. A more refined model, including structural details, has allowed us to verify 
this assumption and will be presented in a future paper || . Transmission Electron Micrographs 
of dilute solutions of side-chain amphiphilic polymers, with and without surfactant, show 
self-assembled hydrophobic microdomains forming locally along the polymeric chains, while 
the chains seem to maintain overall extended conformations [Q] . This observation is explained 
by the highly hydrophilic polymer backbone (usually a polyelectrolyte) and the high density 
of hydrophobic side groups along the chain. We are allowed, therefore, to consider only local 
association along the chain sequence, which greatly simplifies the model. The current work is 
restricted to the dilute limit, where inter-chain effects can be neglected and the surfactant is 
below its critical micelle concentration, in accordance with the relevant experiments. 

Consider a flexible chain of N segments in contact with a surfactant reservoir of chemical 
potential fi s . To each segment we assign a vector, u„, and denote the number of surfactants 
bound to it by (p n . The energy gained by a segment joining a microdomain is denoted by 
a, and the one gained by a surfactant bound to such a microdomain is 7. The energetic 
parameters, a and 7, represent, respectively, measures of the polymer hydrophobicity and 
polymer-surfactant affinity. (All energies in this Letter are in units of k B T.) The partition 
function for the polymer-surfactant system is 

Z ps = Tr [ Vu n exp(-H ps /k B T) 

{Vn}J 

N N-l ^ ^ N 

H ps /k B T = ^A„«^-^ a 4 m (u„ + i-u„) 2 -/i 5 ^y„ (f) 

71 — 1 71 — 1 71— 1 

where the A n are Lagrange multipliers ensuring the constraints (u^) — 1 ((***) denotes a 
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thermal average over all chain configurations). Some of our results can be more easily derived 
imposing infinitely stiff segment constraints (u^ = 1). Yet, other results, in particular the 
effect of polymer configurations on surfactant binding, cannot be obtained as easily using such 
an approach. We have modelled the association property of the polymer by a local tendency 
of the segments to fold, and the preference of the surfactant to join such clusters by a linear 
coupling to this folding (^] . Not explicitly included in eq. (|l|) , interactions between surfactants 
may arise only indirectly through the coupling to chain configurations. Note that the current 
model contains no structural parameters. The hydrophobic microdomains are described as 
featureless clusters of folded polymer segments attracting free surfactants, as illustrated in 
fig. lb. In a more refined model, steric and other repulsive effects should compete with the 
association, resulting in a finite radius of curvature for the microdomains S . Such microscopic 
details, however, do not alter the thermodynamic behaviour which concerns us in this work. 

In the surfactant-free case, the partition function (Q) becomes similar to the one for semi- 
flexible polymers ]Tc| ], 

- N N-l 

Z p = Vu n cxp[- Kul + J (Un+l - u„) 2 ] (2) 

J ra=l n=l 

Yet, because of the tendency of segments in our model to anii-align (a > 0), we need first to 
integrate over the "sub-lattice" U2, 114, . . . before taking the continuum limit. In addition, we 
replace all X n with a single A. (This approximation amounts to a non-extensive correction in 
the free energy, which becomes negligible for N — ► 00 Q.) The resulting expression is 

/ \ 3JV/4 r. ~ p N /j \2 T » TV 

Z p = / V ^ n ) ex P[~? / ( ~T~ ) — u\n)An] (3) 



A - a/2 J J ' 2 J \dn J 2 J 

where A = A(A— a)/(A— a/2), a = q 2 /4(A— a/2), and n : — » N is now a continuous parameter 
along the chain. (Note that the continuum limit implies a large a and, hence, a large a.) The 
path integral in eq. (JsJ) is readily evaluated using an analogy with a three-dimensional quantum 
oscillator and the ground-state dominance p2|, yielding 



3AT/4 

7T N 



Z » = [\^/2 ] 6XP 



-(3N/2)J\/c 



(4) 



Imposing now the constraint on (w 2 ), — d\ogZ p /d\ = N, we find X — a + 0(1 /a), X = 
9/2a+0(l/a 2 ), anda = a/2+0(l/a). The correlation between alternating segment directions 
along the chain is found from eq. (||) as 

(u(n) -u(n')> = exp(-|n-n'|/0) C = a /3 (5) 

The value of £ defines a chain length over which alternating segment directions, u„, u„ + 2 • ■ ■, 
are correlated. Since the other segments within this length, u„ + i,u n+ 3 . . ., are assumed, on 
average, to be anti-aligned with those vectors, we identify £ as the average number of segments 
in a single folded cluster (see fig. lb), i.e., in a single hydrophobic microdomain. Thus, only 
when a is significantly large does the polymer tend to form microdomains by itself. The free 
energy of the chain is 

F p (a) = - hxZ p - XN = N[-a + (3/4) log(a/27r)] (6) 

and can be viewed as composed of competing curvature energy (decreasing with a) and entropy 
(increasing with a) terms. 
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Fig. 2. - (a) Surfactant-polymer binding isotherms. Non-cooperative binding is shown for 
e = 0.1 (solid) and 0.5 (dashed) using eq. (||); cooperative binding is presented for e = 3 
(dotted) using eq. jl^). (b) Diagram presenting the surfactant-polymer binding behaviour. 
The dashed line separates two binding regimes: for small e the binding is non-cooperative, 
whereas for large e the binding is cooperative and exhibits a sharp (though continuous) increase 
at the cac. The values of /x 8 + 7 corresponding to the cac are represented by the solid line. 



We now return to the full polymer-surfactant partition function ([!]). For simplicity let us 
assume that at most one surfactant can bind to a monomer, i.e., tp n is either or 1. Tracing 
over {ip n } and expanding to 2nd order in e e 7/a, we find the following expression for the 
free energy: 



/) + ^-^/(3 + /) £ 2 



(7) 



where / = [1 + c-^+t)]- 1 is the average fraction of bound surfactant if it had been simply 
attracted to the polymer without coupling to its curvature. The first term in eq. (Q) is the 
free energy of the surfactant-free chain, eq. (||), the second describes a regular lattice-gas 
contribution, and the remaining terms represent the surfactant-polymer interaction. From 
eq. (Q) a binding isotherm can be calculated, 



<p ~ f[l - (3/4)(l - f)e + (3/32)(l - /)(3 + 2/)e 2 



(8) 



where (p is the average fraction of bound surfactant. Figure 2a shows two such isotherms (solid 
and dashed lines) for two different, small values of e. In this regime the binding increases 
gradually with fj, s (non-cooperative binding). For larger values of e the isotherms are shifted 
to larger fi s and become steeper. Let us define a chemical potential, /j,*, for which the amount 
of bound surfactant is appreciable, i.e., ip — 1/2. A corresponding cooperativity parameter is 
defined as C = dip/d^l^*. From eq. (JsJ) we readily find 



Ms * ~ - 7 +(3/4)e-(3/8> 
C ~ (l/4)[l + (3/16)e 2 ] 



(9) 
(10) 
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The last, positive contributions to the average binding (|J) and cooperativity ( |l0|) indicate an 
effective short-range attraction between bound surfactants induced by chain configurations. 
Once a molecule has bound to the chain, a folded hydrophobic region is formed, favouring 
binding of another molecule to the same region. Effective attraction is a general feature of 
annealed impurities. The important point, however, is that the attraction strength (resulting in 
binding cooperativity) depends on a single parameter, e, denoting the ratio between surfactant- 
polymer affinity and polymer hydrophobicity. 

For larger e, the value of /i* can be estimated as 



Ms = --j(( U «+l -U") 2 } 



3e 

v =i/2 = ~ 7+ 4(l + e/2) 



To obtain similar estimates for the binding isotherm and cooperativity we return to eq. (Q) 
and make a cumulant expansion about the surfactant-free Hamiltonian before tracing over 
{</?„}. Restricting to two-body interactions between nearest bound surfactants, the resulting 
partition function becomes analogous, to a one-dimensional lattice-gas or Ising model with an 
attraction term proportional to e 2 . (EJ) Calculation of the binding isotherm and cooperativity 
is then straightforward and gives 



9[9-l + V(9-ir+2h] + h _ f ^ hs2ge -^ /B (12) 



(g-l) 2 + (g + l)^/(g-iy + 2h + 2h 1 1-f 
C = (l/4)cxp[(3/16) e 2 ] (13) 

For small e, eqs. (|ll])-(|l3|) coincide with eqs. (f)-©. Equation © is an important result 
since it demonstrates how sensitively the binding cooperativity depends on e. In the regime 
of large e the binding isotherm has a sharp (albeit continuous) slope at /i* , as demonstrated 
by the dotted curve in fig. 2a. Hence, /i s = /i* can be associated with </>*, the so-called critical 
aggregation concentration (cac) of surfactant, fi* — ln0*. 

Our findings are summarized in the diagram drawn in fig. 2b. Depending on the value 
of e, two distinct binding regimes are predicted. The non-cooperative binding regime of 
small e (a 3> 7) corresponds to a situation where the surfactant almost does not affect the 
self-assembly of the polymeric side groups. The polymer is hydrophobic enough by itself 
to form microdomains and the bound surfactants merely join the already existing micelles, 
swelling them a little. On the other hand, the cooperative binding regime of large f (7 » a) 
corresponds to the opposite situation where the surfactant triggers the self-assembly. Above 
the cac, the number of repeat units per hydrophobic microdomain jumps from roughly a to 
about a + 7 ^> a. 

In recent experiments |3| both polymer hydrophobicity and surfactant-polymer affinity were 
changed by controlling the polymer charge and salt concentration (probably affecting both a 
and 7) and the length of the hydrophobic side chains (mainly affecting a) . In those experiments 
the binding cooperativity could, indeed, be sensitively tuned by changing either the polymer 
hydrophobicity or its affinity to the surfactant, and the two distinct binding regimes described 
by our model were clearly observed. 

Similar to regular surfactant micellisation, our self-assembly model yields a well-defined 
surfactant aggregation number determined by the correlation length of the polymer chain (the 
number of repeat units per microdomain) and the average fraction of bound surfactant. In our 
case of surfactant-polymer self-assembly, the aggregate size is restricted by the entropy of the 

(^For large e, many-body interactions among surfactants may become significant. Yet, since in 
our case all such contributions are attractive, they would merely enhance the binding cooperativity 
and further support our findings. 
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amphiphilic backbone, which imposes a finite correlation length along the chain. Moreover, 
like regular miccllisation, the self-assembly at the cue, though sharp, is not a first-order 
transition. In our model this is due to the one-dimensionality of the effective interactions 
between surfactants. 

We have presented a model which accounts for the experimentally observed diverse be- 
haviour of dilute solutions of amphiphilic side-chain polymers and surfactants. Our results 
demonstrate the delicate balance in these mixed systems. By changing the parameters of 
the polymer and surfactant it is possible to cause aggregation and extensively modify the 
macroscopic properties of the solution. While capturing the thermodynamics of self-assembly 
in the system and accounting for the polymer-induced interactions between surfactants, our 
model still lacks the microscopic structural details required for a more complete descrip- 
tion of the assembled microdomains themselves. Such a refinement will be presented in a 
future paper ||. Another important extension of the model is to consider inter-chain as- 
sociation in more concentrated solutions and relate structural properties to rheological ones. 
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